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The problem of flame propagation is studied as an example 
of unstable fronts that wrinkle on many scales. The analytic 
tool of pole expansion in the complex plane is employed to 
address the interaction of the unstable growth process with 
random initial conditions and perturbations. We argue that 
the effect of random noise is immense and that it can never be 
neglected in sufficiently large systems. We present simulations 
that lead to scaling laws for the velocity and acceleration of 
the front as a function of the system size and the level of noise, 
and analytic arguments that explain these results in terms of 
the noisy pole dynamics. 

PACS numbers 47.27.Gs, 47.27.Jv, 05.40. +j 



I. INTRODUCTION 

In our last paper about flame front propagation we left 
two open problems. The first one is how to explain the 
existence of small dependence parameters of problem on 
the noise. The second problem is how to calculate nu- 
merically such values as excess number of poles in system 
, number of poles that appear in the system in unit of 
time, life time of pole. To solve these problem we write 
this paper. 



II. EQUATIONS OF MOTION AND 
POLE-DECOMPOSITION IN THE CHANNEL 
GEOMETRY 

It is known that planar flames freely propagating 
through initially motionless homogeneous combustible 
mixtures are intrinsically unstable. It was reported that 
such flames develop characteristic structures which in- 
clude cusps, and that under usual experimental condi- 
tions the flame front accelerates as time goes on. A model 
in 1 + 1 dimensions that pertains to the propagation of 
flame fronts in channels of width L was proposed in [4]. 
It is written in terms of position h(x, t) of the flame front 
above the rc-axis. After appropriate rescalings it takes 
the form: 



dh(x,t) 
dt 



dh(x, t) 
dx 



d 2 h(x,t) 
dx 2 



I{h(x,t)} + 1 . 

(1) 



The domain is < x < L, v is a parameter and we use 
periodic boundary conditions. The functional I[h(x, t)\ is 
the Hilbert transform of derivative which is conveniently 
defined in terms of the spatial Fourier transform 



h(x,t) = / e lkx h(k,t)dk 
I[h(k,t)] = \k\h(k,t) 



(2) 
(3) 



For the purpose of introducing the pole-decomposition it 
is convenient to rescale the domain to < 8 < 2tt. Per- 
forming this rescaling and denoting the resulting quanti- 
ties with the same notation we have 



dh(0,t) 
dt 



1 

2L 2 



dh(9,t) 



00 



v d 2 h{0,t) 
L 2 dO 2 



L 



I{h{0,t)} + 1 



(4) 



In this equation L = L/2ir. Next we change variables to 
u(9,t) = dh(0,t)/dO. We find 



du{6,t) u(6,t)du(6,t) v d 2 u{0,t) 



dt 



L 2 



00 



L 2 dO 2 



L 



l{u(0,t)} . 
(5) 



It is well known that the flat front solution of this equa- 
tion is linearly unstable. The linear spectrum in k- 
representation is 



u k = \k\/L - vk 2 /L 2 



(6) 



There exists a typical scale k max which is the last unsta- 
ble mode 



k 



(7) 



Nonlinear effects stabilize a new steady-state which is 
discussed next. 

The outstanding feature of the solutions of this equa- 
tion is the appearance of cusp-like structures in the de- 
veloping fronts. Therefore a representation in terms of 
Fourier modes is very inefficient. Rather, it appears very 
worthwhile to represent such solutions in terms of sums of 
functions of poles in the complex plane. It will be shown 
below that the position of the cusp along the front is de- 
termined by the real coordinate of the pole, whereas the 



1 



height of the cusp is in correspondence with the imagi- 
nary coordinate. Moreover, it will be seen that the dy- 
namics of the developing front can be usefully described 
in terms of the dynamics of the poles. Following [8,9,11,7] 
we expand the solutions u(9, t) in functions that depend 
on N poles whose position Zj (t) = Xj (t) + iyj (t) in the 
complex plane is time dependent: 



N 



i(o,t) = cot 

N 



2 S m[e- Xj (t)} 



cosh[t/j(i)] — cos[6 — xj(t)] ' 



(8) 



N 

h{6,t) =2i/^ln[cosh(j/j-(i)) - cos(0 - Xj(t)) +C(t) . 



3=1 



(9) 



In (9) C(t) is a function of time. The function (9) is a 
superposition of quasi-cusps (i.e. cusps that arc rounded 
at the tip). The real part of the pole position (i.e. x 3 ) 
is the coordinate (in the domain [0, 2ir\) of the maximum 
of the quasi-cusp, and the imaginary part of the pole 
position (i.e yj) is related to the depth of the quasi-cusp. 
As yj decreases the depth of the cusp increases. As yj — > 
the depth diverges to infinity. Conversely, when yj —> 
oo the depth decreases to zero. 

The main advantage of this representation is that the 
propagation and wrinkling of the front can be described 
via the dynamics of the poles. Substituting (8) in (5) we 
derive the following ordinary differential equations for the 
positions of the poles: 



dt 



2N , _ 

[u cot {- 

k=i,kjtj 



L 



sign[Im(z 



(10) 



We note that in (8), due to the complex conjugation, we 
have 2N poles which are arranged in pairs such that for 
j < N Zj +N — Zj. In the second sum in (8) each pair of 
poles contributed one term. In Eq.(10) we again employ 
2N poles since all of them interact. We can write the pole 
dynamics in terms of the real and imaginary parts Xj and 
yj . Because of the arrangement in pairs it is sufficient to 
write the equation for either yj > or for yj < 0. Wc 
opt for the first. The equations for the positions of the 
poles read 



^2 dxj 
dt 



N 

v sin (zj 

k=l,k=£j 



X k ) 



[cosh(yj - y k ) (11) 



-cos( Xj -x k )] 1 + [cosh(?/j + y k ) - cos(xj - x k )] 1 



-2 d V] 



N 



sinh(2/j - y k ) 



lU ' kJtk^j ^ cosh (% - yk) - cos^ - x k ) 



smh(yj + y k ) 



cosh(y :) - + yk) - cos( Xj - x k ) 



+ ucoth(yj) — L. (12) 



We note that if the initial conditions of the differential 
equation (5) are expandable in a finite number of poles, 
these equations of motion preserve this number as a func- 
tion of time. On the other hand, this may be an unstable 
situation for the partial differential equation, and noise 
can change the number of poles. This issue will be exam- 
ined at length in Section III. We turn now to a discussion 
of the steady state solution of the equations of the pole- 
dynamics. 



A. Qualitative properties of the stationary solution 

The steady-state solution of the flame front propagat- 
ing in channels of width 2ir was presented in Ref. [9]. 
Using these results we can immediately translate the dis- 
cussion to a channel of width L. The main results are 
summarized as follows: 

1 . There is only one stable stationary solution which is 
geometrically represented by a giant cusp (or equiv- 
alently one finger) and analytically by N(L) poles 
which are aligned on one line parallel to the imag- 
inary axis. The existence of this solution is made 
clearer with the following remarks. 

2. There exists an attraction between the poles along 
the real line. This is obvious from Eq.(ll) in 
which the sign of dxj/dt is always determined by 
sm(xj — Xk). The resulting dynamics merges all the 
x positions of poles whose y-position remains finite. 

3. The y positions are distinct, and the poles are 
aligned above each others in positions yj-i < yj < 
yj+i with the maximal being yjv(L)- This can be 
understood from Eq.(12) in which the interaction 
is seen to be repulsive at short ranges, but changes 
sign at longer ranges. 

4. If one adds an additional pole to such a solution, 
this pole (or another) will be pushed to infinity 
along the imaginary axis. If the system has less 
than N(L) poles it is unstable to the addition of 
poles, and any noise will drive the system towards 
this unique state. The number N(L) is 



JV(L) 



-(- 

2\v 



1 



(13) 



where 



is the integer part. To see this con- 



sider a system with N poles and such that all the 
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values of yj satisfy the condition < yj < ymax- 
Add now one additional pole whose coordinates are 
z a = [x a ,y a ) with y a > y ma x- From the equa- 
tion of motion for y a , (12) we see that the terms in 
the sum are all of the order of unity as is also the 
cot(y a ) term. Thus the equation of motion of y a is 
approximately 



dy a 
dt 



2N+1 
L 2 



1 

L 



(14) 



The fate of this pole depends on the number of 
other poles. If N is too large the pole will run 
to infinity, whereas if N is small the pole will be 
attracted towards the real axis. The condition for 
moving away to infinity is that N > N(L) where 
N(L) is given by (13). On the other hand the y 
coordinate of the poles cannot hit zero. Zero is 
a repulsive line, and poles are pushed away from 
zero with infinite velocity. To see this consider a 
pole whose yj approaches zero. For any finite L 
the term coth(j/j) grows unboundedly whereas all 
the other terms in Eq.(12) remain bounded. 

5. The height of the cusp is proportional to L. The 
distribution of positions of the poles along the line 
of constant x was worked out in [9] . 

We will refer to the solution with all these properties as 
the Thual-Frisch-Henon (TFH)-cusp solution. 
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FIG. 1. The 

dependence of the average velocity v on the system size L for 
/°- 5 = 0, 2.7e-6, 2.7e-5, 2.7e-4, 2.7e-3, 2.7e-2, 2.7e-l, 0.5, 1.3, 2.7. 

and the Fourier amplitudes rj k are correlated according 
to 



<r ]k (t)r ] * k ,(t') >=£-5 kik ,6(t-t') 



(16) 



We will first examine the result of numerical simulations 
of noise-driven dynamics, and later return to the theo- 
retical analysis. 



III. ACCELERATION OF THE FLAME FRONT, 
POLE DYNAMICS AND NOISE 

A major motivation of this Section is the observation 
that in radial geometry the same equation of motion 
shows an acceleration of the flame front. The aim of 
this section is to argue that this phenomenon is caused 
by the noisy generation of new poles. Moreover, it is our 
contention that a great deal can be learned about the 
acceleration in radial geometry by considering the effect 
of noise in channel growth. In Ref. [9] it was shown that 
any initial condition which is represented in poles goes to 
a unique stationary state which is the giant cusp which 
propagates with a constant velocity v = 1/2 up to small 
1/L corrections. In light of our discussion of the last sec- 
tion we expect that any smooth enough initial condition 
will go to the same stationary state. Thus if there is no 
noise in the dynamics of a finite channel, no acceleration 
of the flame front is possible. What happens if we add 
noise to the system? 

For concreteness we introduce an additive white-noise 
term rj{9, t)to the equation of motion (5) where 

V (6,t) =£»fe(i)exp(ifc0) , (15) 

fe 



A. Noisy Simulations 

Previous numerical investigations [6,15] did not intro- 
duce noise in a controlled fashion. We will argue later 
that some of the phenomena encountered in these simu- 
lations can be ascribed to the (uncontrolled) numerical 
noise. We performed numerical simulations of Eq.(5 us- 
ing a pseudo-spectral method. The time-stepping scheme 
was chosen as Adams-Bashforth with 2nd order presi- 
cion in time. The additive white noise was generated in 
Fourier-space by choosing r) k for every k from a flat distri- 
bution in the interval [— \J^^- \J 2^-] • We examined the 
average steady state velocity of the front as a function of 
L for fixed / and as a function of / for fixed L. We found 
the interesting phenomena that are summarized here: 

1. In Fig. 2 we can see two different regimes of behavior 
the average velocity v as function of noise /°.5 for 
fixed system size L. For the noise / smaller then 
same fixed value f cr 

v~ft. (17) 

For these values of / this dependence is very weak, 
and £ w 0.02. For large values of / the dependence 
is much stronger 
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FIG. 2. The dependence of the average velocity v on the 
noise f 5 for L=10, 40, 80. 

2. In Fig.l we can see growth of the average velocity 
v as function of the system size L. After some val- 
ues of L we can see saturation of the velocity. For 
regime / < f cr the growth of the velocity can be 
written as 

u~L", n « 0.40 ± 0.05 . (18) 




FIG. 3. The dependence of the cusps positions on 
time.L = 80 v = 0.1 / = 9e - 6 

(ii) Regime II Strong dependence of poles number 
SN on noise /. 

(iii) Regime III Saturation poles number SN on 
noise /, so we see very small dependence of this 
number on noise 

SN ~ J 003 (19) 



B. Calculation of Poles Number in the System 

The interesting problem that we would like to solve 
here it is to find number of poles that exist in our system 
outside the giant cusp. We can make it by next way: to 
calculate number of cusps (points of minimum or inflex- 
ional points) and their position on the interval 6 : [0, 2-7r] 
in every moment of time and to draw positions of cusp 
like function of time, see Fig. 3. 

We assume that our system is almost all time in " quasi- 
stable" state, i.e. every new cusp that appears in the 
system includes only one pole. By help pictures obtained 
by such way we can hnd 

1. By calculation number of cusp in some moment of 
time and by investigation of history of every cusp 
(except the giant cusp) , i.e. how many initial cusps 
take part in formation this cusp, after averaging 
with respect to different moments of time we can 
find mean number of poles that exist in our system 
outside the giant cusp. Let us denote this number 
SN. We can see four regimes that can be define 
with respect to dependence of this number on noise 
/: 

(i) Regime I: Such small noise that no poles exist 
in our system outside the giant cusp. 



Saturated value of SN is defined by next formula 
*JV«tf(Z)/2«i£ (20) 



where N(L) « |~ is number of poles in giant cusp. 

(iv) Regime IV We again see strong dependence of 
poles number SN on noise /. 



SN 



o.i 



(21) 



Because of numerical noise we can see in most of 
simulations only regime III and IV. In future if we 
don't note something different we discuss regime 
III. 

2. By calculation of new cusp number we can find 
number of poles that appear in the system in unit 
of time Of. In regime III 



dN 
~dt 



f0.03 



(22) 



Dependence on L and v define by 
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FIG. 5. The dependence of excess poles number 8N on the FIG. 7. The dependence of excess poles number SN on the 

noise f 5 . v = 0.1 L = 40, 80. system size L. v = 0.1 f 5 = 9e - 6 



(23) 
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FIG. 8. The dependence of poles number in time unit 
dN/dt on the parameter v. L = 80 v = 0.1. 
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FIG. 9. The dependence of excess poles number SN on the 
the parameter v. L — 80 v — 0.1. 



dN _1_ 

It^^ 2 



And in regime IV 



dN 



fO.l 



(24) 



(25) 



C. Theoretical Discussion of the Effect of Noise 

1. The Threshold of Instability to Added Noise. Transition 
from regime I to regime II 

First we present the theoretical arguments that explain 
the sensitivity of the giant cusp solution to the effect of 
added noise. This sensitivity increases dramatically with 
increasing the system size L. To see this we use again 
the relationship between the linear stability analysis and 
the pole dynamics. 

Our additive noise introduces perturbations with all 
fc-vectors. We showed previously that the most unsta- 
ble mode is the k = 1 component ^4isin(#). Thus the 
most effective noisy perturbation is r]i sin(0) which can 
potentially lead to a growth of the most unstable mode. 
Whether or not this mode will grow depends on the am- 
plitude of the noise. To see this clearly we return to the 
pole description. For small values of the amplitude A\ 
we represent A\ sin(#) as a single pole solution of the 
functional form ve~ y sine?. The y position is determined 
from y = — log|Ai|/V, and the ^-position is 9 = n for 
positive A\ and 9 = for negative A\. From the analysis 
of Section III we know that for very small A\ the fate 
of the pole is to be pushed to infinity, independently of 
its 9 position; the dynamics is symmetric in A\ — > — A\ 
when y is large enough. On the other hand when the 
value of Ai increases the symmetry is broken and the 9 
position and the sign of A\ become very important. If 
A\ > there is a threshold value of y below which the 
pole is attracted down. On the other hand if A\ < 0, 
and 9 = the repulsion from the poles of the giant cusp 
grows with decreasing y. We thus understand that qual- 
itatively speaking the dynamics of Ai is characterized by 
an asymmetric "potential" according to 



A, =- 



dA 1 ' 
V{Ax) = XA\-aA\ + .. 



(26) 
(27) 



^From the linear stability analysis we know that A w 
v/L 2 , cf. Eq.(14). We know further that the threshold 
for nonlinear instability is at A\ w v z /L 2 , cf. Eq(??). 



G 



This determines that value of the coefficient a « 2/iv 2 . 
The magnitude of the "potential" at the maximum is 



V(A max ) « is 7 /L 6 



(28) 



The effect of the noise on the development of the mode 
Ai sin 9 can be understood from the following stochastic 
equation 



dV(A 1 ) 

dAy 



(29) 



It is well known [17] that for such dynamics the rate of 
escape R over the "potential" barrier for small noise is 
proportional to 



R 



L 2 



exp 



(30) 



The conclusion is that any arbitrarily tiny noise becomes 
effective when the system size increase and when v de- 
creases. If we drive the system with noise of amplitude 
£ the system can always be sensitive to this noise when 
its size exceeds a critical value L c that is determined by 
J- v 7 / L®. This formula defines transition from regime 
I (no poles) to regime II. For L > L c the noise will in- 
troduce new poles into the system. Even numerical noise 
in simulations involving large size systems may have a 
macroscopic influence. 

The appearance of new poles must increase the veloc- 
ity of the front. The velocity is proportional to the mean 
of (u/L) 2 . New poles distort the giant cusp by additional 
smaller cusps on the wings of the giant cusp, increasing 
u 2 . Upon increasing the noise amplitude more and more 
smaller cusps appear in the front, and inevitably the ve- 
locity increases. This phenomenon is discussed quantita- 
tively in Section III. 



2. Verifing of asymmetric "potential" form 

If we know distribution of poles in the giant cusp we 
can find form of "potential" and verify expressions for 
A, A max and 9V q^ ■ For measurement of A\ we use 

formula A x = Ive'V and = 4i/|jfe-» where |f 

can be find from equation motion of poles. Numerical 
measurements were made for L — 2nv,n — integer, n > 2. 
For graphs that we use for fixed v and variable L\v = 
0.005 L : [1,150]. For fixed L and variable v; L = \ 
and v : [0.005,0.05]. We can find A max as zero-point of 
9 . Obtain results are next. For Ama % L as function 
of L (Fig 10) 
and v (Fig 11) 

we can see that — is almost constant. Depen- 
dence ^" (A/v(L) i s defined by position of the upper 
pole.)as function L (Fig 12) 

and v (Fig 13) 
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FIG. 10. Dependence of the normalized amplitude 
A max L 2 /v 3 on the system size L. 
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FIG. 11. Dependence of the normalized amplitude 
AmaxL 2 / v 3 on the parameter v. 
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FIG. 12. Relationship between amplitude defined by the 
minimum of potential A max and amplitude defined by posi- FI G. 14. Dependence of the normalized parameter XL 2 /i> 

tion of upper pole A N ( L ) as function of the system size L. on tne system size L. 
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FIG. 13. Relationship between amplituda defined by the 
minimum of potential A max and amplituda defined by posi- 
tion of upper pole A N ^ L ) as function of the the parameter 
v. 
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FIG. 15. Dependence of the normalized parameter XL 2 jv 
on the parameter v. 
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FIG. 16. The first odd eigenfunction obtained from tradi- 
tional stability analysis. 



is almost Const « 0.5. From graphs Fig 14 (for L) 
and Fig 15 

(for v) demonstrate us that ^j- is almost constant. 
Fig 16 give us full dependence of £ on L. These results 
are in good agreement with the theory. 



3. The Noisy Steady State and its Collapse with Large Noise 
and System Size 

In this subsection we discuss the response of the giant 
cusp solution to noise levels that are able to introduce a 
large number of excess poles in addition to those exist- 
ing in the giant cusp. We will denote the excess number 
of poles by SN. The first question that we address is 
how difficult is it to insert yet an additional pole when 
there is already a given excess SN. To this aim we esti- 
mate the effective potential Vsn{A\) which is similar to 
(27) but is taking into account the existence of an excess 
number of poles. A basic approximation that we employ 
is that the fundamental form of the giant cusp solution 
is not seriously modified by the existence of an excess 
number of poles. Of course this approximation breaks 
down quantitatively already with one excess pole. Qual- 
itatively however it holds well until the excess number of 
poles is of the order of the original number N(L) of the 
giant cusp solution. Another approximation is that the 
rest of the linear modes play no role in this case. At this 
point we limit the discussion therefore to the situation 
SN < N(L) (regime II). 

To estimate the parameter A in the effective potential 
we consider the dynamics of one pole whose y position 
y a is far above y ma x- According to Eq.(14) the dynamics 
reads 



dy a _ 2v{N{L) +5N) _ 1 
dt L 2 L 



(31) 



Since the N(L) term cancels against the L _1 term (cf. 
Sec. II A), we remain with a repulsive term that in the 
effective potential translates to 



A = 



vSN 



(32) 



Next we estimate the value of the potential at the break- 
even point between attraction and repulsion. In the last 
subsection we saw that a foreign pole has to be inserted 
below y ma x in order to be attracted towards the real axis. 
Now we need to push the new pole below the position 
of the existing pole whose index is N(L) — SN. This 
position is estimated as in Sec III C by employing the 
TFH distribution function (??). We find 



usn 



2 In 



AL 



2 v8N 



(33) 



As before, this implies a threshold value of the amplitude 
of single pole solution A max sin 9 which is obtained from 
equating A max = ve~ VSN . We thus find in the present 
case A max ~ v 3 (5N) 2 /L 2 . Using again a cubic represen- 
tation for the effective potential we find a = 2/(3v 2 SN) 
and 



V(A ri 



_ 1 v 7 {SNf 
c) ~ 3 U> 



(34) 



Repeating the calculation of the escape rate over the po- 
tential barrier we find in the present case 



R 



;SN 



L 2 



■ cxp 



-S(8N)>/iL' 



(35) 



For a given noise amplitude £ there is always a value 
of L and v for which the escape rate is of O(l) as long 
as SN is not too large. When SN increases the escape 
rate decreases, and eventually no additional poles can 
creep into the system. The typical number SN for fixed 
values of the parameters is estimated from equating the 
argument in the exponent to unity: 



SN 



L e /v 7 



1/5 



(36) 



We can see that SN depend on noise / very seriously. 
It is not the case in regime III. Let us find conditions 
of transition from regime II to regime III, where we see 
saturation of SN with respect to noise /. 

(i) We use for the amplitude of pole solution that really 



equal to 



2u sin 9 



expression A„ 



Ave~ VSN but 



cosh(y5 jy ) — COS 9 

it is right only for big number ysN- For ysN < 1 better 
approximation is A max = -4r-- From equation (33) we 
can find that the boundary value ysN — 1 correspond to 
SN « N(L)/2 
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(ii) We use expression ysN ~ 2 In 



4L 



8eL 

t 2 v{N(L)-8N) 



^U5N\ but for bi § 
value of SN better approximation that can be find the 

same way is ysN ~ ! ^t{N(L) — 8N)ln 

These expressions give us nearly equal result for SN 
N(L)/2. 

From (i) and (ii) we can make next conclusions 

(a) Transition from regime II to regime III happens for 
nearly SN w N(L)/2 

(b) Using new expression in (i) and (ii) for amplitude 
A m ax and ysN we can find for noise ^ in regime III: 



noise and instability and which results in the acceleration 
of the flame front. 

To estimate the velocity of the flame front we need 
to write down an equation for the mean of < dh/dt > 
given an arbitrary number N of poles in the system. This 
equation follows directly from (4): 



dh 
dt 



i i r 2 



i 2 d9 



(42) 



After substitution of (8) in (42) we get, using (11) and 
(12) 



/ 
L 



V(A max ) ~ \A\ 



vSN 4^ 2 



,57V 



v (N(L) - SN) 4 



(37) 



This expression define very slow dependence of SN on 
noise £ for SN > N(L)/2 that explain noise saturation 
of SN for regime III. 

(c) Form of the giant cusp solution is defined by poles 
that are closely to zero with respect to y. For regime 
III N(L)/2 poles that have position y < ysN=N(L)/2 = 1 
stay on these place. This result explain why giant cusp 
solution is not seriously modified for regime III. 

From Eq.(36) by help of boundary condition 

5N^N{L)/2 (38) 
boundary noise fb between regime II and III can be found 



fb ~ v 1 



(39) 



The basic equation describing pole dynamics is next 
dN _ SN 



(40) 



where ^ is number of poles that appear in unit of 
time in our system, SN is excess number of poles, T is 
mean life time of pole (between appearing and merging 
with giant cusp). Using result of numerical simulations 
for ^ and (38) we can find for T 



T 



SN 
In 

dt 



0.2 



(41) 



So life time proportional to v and depend on system 
size L very weakly. 



D. The acceleration of the flame front due to noise 

In this section we estimate the scaling exponents that 
characterize the velocity of the flame front as a function of 
the system size. Our arguments in this section are even 
less solid than the previous ones, but nevertheless we 
believe that we succeed to capture some of the essential 
qualitative physics that underlies the interaction between 



dt V L 



v 2 N< 
~1? 



(43) 



The estimates of the second and third terms in this equa- 
tion are straightforward. Writing TV = N(L)+SN(L) and 
remembering that N(L) ~ Ljv and SN(L) ~ N(L)/2 
we find that these terms contribute 0(1). The first term 
will contribute only when the current of poles is asym- 
metric. Since noise introduces poles at a finite value of 
Vmin, whereas the rejected poles stream towards infin- 
ity and disappear at boundary of nonlinearity defined by 
position of highest pole 

4L 



2 In 



(44) 



, we have an asymmetry that contributes to the veloc- 
ity of the front. To estimate the first term let us define 



, l+dl , 

dy k _ m-^ dy k 
a[ 2^ dt > - 2^ At 



(45) 



where Y^\ +dl is sum over poles that are on the in- 



terval y : [1,1 + dl]. We can write 



dt 



Where d(%2 ~dt)up nux °f P°l cs moving up and 

d{J2 ^§t)down flux of poles moving down. 

For these flux we can write 

, { sr^dy k ,,\^dy k dN 
d(}_^ -jr)u P , d{}^ —)down < -j-dl (47) 



dt 

So for the first term 
N a 

o< y^ = 



fc=i 



(48) 



f 

J Vn 



dl 



I Vmin 

.dN 

dt ^ max ^ n 
dN 

S TrV max 

dt 



dl 



dl 
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Because of slow(ln) dependence of y max on L and v ^ 
term define oder of nonlinearity for first term. This term 
zero for symmetric current of poles and achieves max- 
imum for maximal asymmetric current of poles. Com- 
parison v ~ L°- 42 /°- d and ^ ~ L o.8y0.03 confirm this 
calculation. 

IV. SUMMARY AND CONCLUSIONS 

The main two messages of this paper are: (i) There is 
an important interaction between the instability of de- 
veloping fronts and random noise; (ii) This interaction 
and its implications can be understood qualitatively and 
sometimes quantitatively using the description in terms 
of complex poles. 

The pole description is natural in this context firstly 
because it provides an exact (and effective) representa- 
tion of the steady state without noise. Once one succeeds 
to describe also the perturbations about this steady state 
in terms of poles, one achieves a particularly transparent 
language for the study of the interplay between noise and 
instability. This language also allows us to describe in 
qualitative and semi-quantitative terms the inverse cas- 
cade process of increasing typical lengths when the sys- 
tem relaxes to the steady state from small, random initial 
conditions. 

The main conceptual steps in this paper are as follows: 
firstly one realizes that the steady state solution, which is 
characterized by N{L) poles aligned along the imaginary 
axis is marginally stable against noise in a periodic array 
of L values. For all values of L the steady state is non- 
linearly unstable against noise. The main and foremost 
effect of noise of a given amplitude / is to introduce an 
excess number of poles SN(L, /) into the system. The 
existence of this excess number of poles is responsible for 
the additional wrinkling of the flame front on top of the 
giant cusp, and for the observed acceleration of the flame 
front. By considering the noisy appearance of new poles 
we rationalize the observed scaling laws as a function of 
the noise amplitude and the system size. 

Theoretically we therefore concentrate on estimating 
5N(L, /). The measurements do not test our theoretical 
consideration directly, but rather test the dependence of 
the velocity on L and /. The only direct test for our 
theory is the critical line shown in Fig. 7. The measured 
exponent is in accord with our analytic estimates. Nev- 
ertheless we note that some of our consideration are only 
qualitative. For example, we estimated SN(L, /) by as- 
suming that the giant cusp solution is not seriously per- 
turbed. On the other hand we find a flux of poles going 
to infinity due to the introduction of poles at finite values 
of y by the noise. The existence of poles spread between 
Vmax and infinity is a significant perturbation of the gi- 
ant cusp solution. Thus also the comparison between the 
various scaling exponents measured and predicted must 



be done with caution; we cannot guarantee that those 
cases in which our prediction hit close to the measure- 
ment mean that the theory is quantitative. However we 
believe that our consideration extract the essential ingre- 
dients of a correct theory. 

The "phase diagram" as a function of L and / in this 
system consists of three regimes. In the first one, dis- 
cussed in Section III C 1 , the noise is too small to have 
any effect on the giant cusp solution. In the second the 
noise introduces excess poles that serve to decorate the 
giant cusp with side cusps. In this regime we find scaling 
laws for the velocity as a function of L and / and we are 
reasonably successful in understanding the scaling expo- 
nents. In the third regime the noise is large enough to 
create small scale structures that are not neatly under- 
stood in terms of individual poles. It appears from our 
numerics that in this regime the roughening of the flame 
front gains a contribution from the the small scale struc- 
ture in a way that is reminiscent of stable, noise driven 
growth models like the Kardar-Parisi-Zhang model. 

One of our main motivations in this research was to 
understand the phenomena observed in radial geometry 
with expanding flame fronts. A full analysis of this prob- 
lem cannot be presented here. We note however that 
many of the insights offered above translate immediately 
to that problem. Indeed, in radial geometry the flame 
front accelerates and cusps multiply and form a hierar- 
chic structure as time progresses. Since the radius (and 
the typical scale) increase in this system all the time, new 
poles will be added to the system even by a vanishingly 
small noise. The marginal stability found above holds 
also in this case, and the system will allow the intro- 
duction of excess poles as a result of noise. The results 
discussed in Rcf. [7] can be combined with the present in- 
sights to provide a theory of radial growth. This theory 
will be offered in a forthcoming publication. 

Finally, the success of this approach in the case of flame 
propagation raises hope that Laplacian growth patterns 
may be dealt with using similar ideas. A problem of 
immediate interest is Laplacian growth in channels, in 
which a finger steady-state solution is known to exist. It 
is documented that the stability of such a finger solution 
to noise decreases rapidly with increasing the channel 
width. In addition, it is understood that noise brings 
about additional geometric features on top of the finger. 
There are enough similarities here to indicate that a care- 
ful analysis of the analytic theory may shed as much light 
on that problem as on the present one. 
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